Heat Transfer between Weakly Coupled Systems: Graphene on a-Si02 
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We study the heat transfer between weakly coupled systems with flat interface. We present a 
simple analytical result which can be used to estimate the heat transfer coefficient. As an application 
we consider the heat transfer between graphene and amorphous Si02. The calculated value of the 
heat transfer coefficient is in good agreement with the value deduced from experimental data. 



Almost all surfaces in Nature and Technology have 
roughness on many different length scales [Ij. When two 
macroscopic solids are brought into contact, even if the 
applied force is very small, e.g., just the weight of the 
upper solid block, the pressure in the asperity contact 
regions can be very high, usually close to the yield stress 
of the (plastically) softer solid. As a result good ther- 
mal contact may occur within each microscopic contact 
region, but owing to the small area of real contact the 
(macroscopic) heat transfer coefficient may still be small. 
In fact, recent studies have shown that in the case of sur- 
faces with roughness on many different length scales, the 
heat transfer is independent of the area of real contact [3|- 
We emphasize that this remarkable and counter-intuitive 
result is only valid when roughness occur over several 
decades in length scale. 

For nanoscale systems the situation may be very dif- 
ferent. Often the surfaces are very smooth with typi- 
cally nanometer (or less) roughness on micrometer-sized 
surface areas, and because of adhesion the solids often 
make contact over a large fraction of the nominal con- 
tact area. The heat transfer between solids in perfect 
contact is usually calculated using the so called diffusive 
mismatch model|, where it is assumed that all phonons 
scatter diffusively and elastically at the interface between 
two materials. In this model there is no direct reference 
to the nature of the solid-solid interaction accross the 
interface, and the model cannot describe the heat flow 
between weakly interacting solids. 

Here we will discuss the heat transfer across perfectly 
flat interfaces, when the interaction between the solids 
is very weak, e.g., of the Van der Waals type, as for 
graphene or carbon nanotubes on many substrates. We 
present a simple analytical result which can be used to es- 
timate the heat transfer coefficient. As an application we 
consider the heat transfer between graphene and amor- 
phous Si02. The calculated value of the heat transfer 
coefficient is in good agreement with the value deduced 
from experimental data. 

Consider the interface between two solids, and assume 
that local thermal equilibrium occurs everywhere except 
at the interface. The energy flow (per unit area) through 
the interface is given byl2l| 

J = a(To - Ti), 



solid 




FIG. 1: A membrane (solid 0) in contact with a semi-infinite 
solid 1. The interfacial surface separation is the sum of the 
equilibrium separation deq and the difference in the surface 
displacements uq—ui, due to thermal movements, where both 
uo and ui are positive when the displacement point along the 
z-axis towards the interior of solid 1. 



where Tq and Ti are the local temperatures at the in- 
terface in solid and 1, respectively. We now present 
a calculation of the heat transfer coefficient a under the 
assumption of very weak coupling between the two solids. 
The stress or pressure acting on the surface of solid 1 
from solid can be written as 

a{:>c,t) = K[uo{x,t) ~ ui(x, i)], 

where uq and ui are the (perpendicular) surface displace- 
ment of solid and 1 (see Fig. [1]) , respectively, and where 
K is a spring constant per unit area characterizing the 
interaction between the two solids. If we define 
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we get 

cr(q,tj) ^ K[uo{c[,uj) -■ui(q,a;)]. 
If solid 1 is semi-infinite we havejj] 

ui{q,uj) = A/(q,w)cr(q,tj), 



(1) 



(2) 



where Af(q, a;) is determined by the elastic properties of 
solid 1. Combining (1) and (2) gives 
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The energy transferred to solid 1 from solid during the 
time period to can be written as 

AE = / (Pxdt ?ii(x,t)cr(x,i), 

where ii — du/dt. One can also write 

AE ^ {2iif I (fqduj (-ia;)ui(q,a;)cr(-q, -w). (4) 

Using (3) and (4) we obtain 



AE = (27r)3 f (fqdiu ^ 



u:KHmM{ci,uj) 



+ KM{ci,Lo)\ 



(|uo(q,w)|^), 



where we have performed an ensemble (or thermal) av- 
erage denoted by (..}. Next, note that 



/I I m2\ ^0^0 

(K(q,^)| ) = J^ 

where Aq is the surface area, and 
1 
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is the displacement correlation function. We can also 
write 



C„„(q,a;) = / d^q'duj' ('Uo(-q, -w)Mo(q', w'))- (7) 

Substituting (6) in (5) gives the heat current Jo^i = 
AE/Ai^ti^ from solid to soHd 1: 

Jo^i =2 d^q dw — — — — ^^i^C'™(q,c^;To). 



|l + i^M(q,Lj)|2 

(8) 
where Cuu depends on the temperature Tq of solid 0. A 
similar equation with Tq replaced by Ti gives the energy 
transfer from solid 1 to solid 0, and the net energy flow 

J = Jo^l — "/l-s-O- 

With graphene in mind, we now assume that solid is 
a membrane or two-dimensional (2D) system. We assume 
that uo(x, t) satisfies 

P^-Q^ = -'^^ V uo - pTi— + /, (9) 

where k w 1 eV [5[ is the bending elasticity, 77 a phe- 
nomcnological friction coefficient, and /(x, i) a stochas- 
tic fluctuating force related to the temperature Tq and 
the friction rj via the fluctuation-dissipation theorem: 

(/(q,^)/(q', ^')> = 2[2T,)--'pi^knTa5[uj + c^')'5(q + q')- 

(10) 
where po = nQm^ is the mass density per unit area of 
the 2D-system (mo is the atom mass and no the number 
of atoms per unit area). On the right hand side of (9) 



there should in principle be another term i^V^Mo, due 
to the interaction with the substrate wall, but because 
of the assumed small magnitude of K, this term can be 
neglected for the relevant wavevectors q. The formalism 
above assumes that the system can be treated classically 
which is the case only if ?iw << fceTb. If this condition is 
not satesfied one can take into account the most impor- 
tant quantum mechanical effects by replacing fee To by 
n(w) == huj [exp(;ia;/fcBTo) - Vf^'"^ . 
From (9) we get 



uo(q,w) 



/(q.'^) 



-puj^ + nq"^ — prjiuj 



(5) Using (7) and (10) we get 



C„„(q,cj) 



2{2TT)-^pijkBTo 
Ipuj'^ — nq'^ + prjiuj\ 



(6) In the limit 77 — > we get 
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where wi = {k/ po)^''^q'^ = c{q)q, where we have defined 
the velocity c{q) = (n/poy^^q. 

Substituting (11) in (8) and assuming weak coupling 
between the solids (i.e., K is small), we get 



Jo- 



27r/9o 



nOO 

/ dq —KHmM{q,uJi). 
Jo ^1 



The heat transfer coefficient a — {Jo-fi — Ji->-o)/{To — Ti) 
is given by 



^ / dq -^KHmM{q,u:i) 
2ttpo Jo wi 



(12) 



Using the expression for Af(q, cj) derived in \4, \&\ and 
wi = c{q)q gives 



where 
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where cl, ct and p\ are the longitudinal and transverse 
sound velocities, and the mass density, respectively, of 
solid 1. The cut off wavevector ^c ~ 7r/ai (a\ is the lat- 
tice constant, or the average distance between two nearby 
atoms) of solid 1 . 



There are two contributions to the integral £,. One is 
derived from c{q) > ci, , but for graphene on a-Si02 stud- 
ied below, this gives only ~ 10% of the contribution to 
the integral. For c(q) < c\^ the term after the Re operator 
is purely imaginary (and will therefore not contribute to 
the integral), except for the case where the denominator 
vanish. It is found that this pole-contribution gives the 
main contribution {^ 90%) to the integral, and corre- 
sponds to the excitation of a Rayleigh surface (acoustic) 
phonon of solid 1. This process involves energy exchange 
between a bending vibrational mode of the graphene and 
a Rayleigh surface phonon mode of solid 1. The denom- 
inator vanish when c{q) = cr where 
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= 0. 



Note that the Rayleigh velocity cr < ct but close to ct. 
For example, when cl/ct = 2, cr, « 0.93ct, and the pole 
contribution to the integral in ^ is 0.083. 

In the model above the heat transfer between the solids 
involves a single bending mode of the membrane or 2D- 
system. In reality there will always be some roughness 
at the interface which will somewhat blurre the wavevec- 
tor conservation rule. We therefore expect some narrow 
band of bending modes to be involved in the energy trans- 
fer, rather than a single mode. Nevertheless, the model 
study above assumes implicitly that, due to lattice non- 
linearity (and defects) , there exist phonon scattering pro- 
cesses which rapidly transfer energy to the bending mode 
involved in the heat exchange with the substrate. This 
requires very weak coupling to the substrate, so that the 
energy transfer to the substrate is so slow that the bend- 
ing mode can be re-populated by phonon scattering pro- 
cesses in the 2D-system, e.g., from the in-plane phonon 
modes, in such a way that its population is always close to 
what would be the case if complete thermal equilibrium 
occurs in the 2D-system. This may require high temper- 
ature in order for multi-phonon scattering processes to 
occur by high enough rates. 

Graphene, the recently isolated 2D-carbon mate- 
rial with unique properties due to its linear electronic 
dispersion, is being actively explored for electronic 
applications [7|. Important properties are the high mobil- 
ities reported especially in suspended graphene, the fact 
that graphene is the ultimately thin material, the stabil- 
ity of the carbon-carbon bond in graphene, the ability to 
induce a bandgap by electron confinement in graphene 
nanoribbons, and its planar nature, which allows estab- 
lished pattering and etching techniques to be applied. 
Recently it has been found that the heat generation in 
graphene field-effect transistors can result in high tem- 
perature and device failure (Sj. Thus, it is important to 
understand the the mechanisms which influence the heat 
flow. 

The graphene — a-Si02 interaction is probably of the 
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FIG. 2: The calculated graphene — a-Si02 interaction energy 
U{d) per graphene carbon atom, as a function of the separa- 
tion d (in nm) between the center of a graphene carbon atom 
and the center of the first layer of substrate atoms. See text 
for details. 



van Der Waals type. In Ref. [9| the interaction be- 
tween the graphene C-atoms and the substrate Si and O 
atoms was assumed to be described by the Lennard- Jones 
(LJ) pair-potentials with different parameters. Here we 
use a simplified picture where the substrate atoms form 
a simple cubic lattice with the lattice constant deter- 
mined by oi = (m/pi)^''^ sa 0.25 nm, where to = 
(TOsi -I- 2mo)/3 ~ 3.32 x 10^^^ kg is the average sub- 
strate atomic mass, and pi « 2200 kg/m'^ the mass den- 
sity of a-Si02. We also use the effective LJ energy pa- 
rameter, e = (esi + 2eo)/3 ~ 5.3 meV, and the bond- 
length parameter a = (csi + 2cro)/3 ~ 0.31 nm. With 
these parameters we can calculate the graphene— a-Si02 
interaction energy, U{d), per graphene carbon atom, as a 
function of the separation d (in nm) between the center of 
a graphene carbon atom and the center of the first layer 
of substrate atoms. We find (see Fig. [5]) the graphene— a- 
Si02 binding energy E'b — —U{dcq) ~ 35 meV per car- 
bon atom, and the force constant k = U"{dcq) (where 



J-cq 



0.32 nm is the equilibrium separation) k 



Kal 



2.4 N/m per carbon atom. This gives the perpendic- 
ular graphene — a-Si02 (uniform) vibration frequency 
cjj^ « [k/moY''^ « 55 cm^"'^, which is similar to what is 
observed for the perpendicular vibrations of linear alkane 
molecules on many surfaces (e.g., about 55 — 60 cm^^ for 
alkanes on metals and on hydrogen terminated diamond 
C(lll)[i3). Using K = k/al = 1.82 x 10^° N/m^, and 
the transverse and longitudinal sound velocities of solid 
1 (ct = 3743 m/s and cl = 5953 m/s), from (13) we 
obtain a « 3 x 10^ W/Km^. 

The heat transfer coefficient between graphene and a 
perfectly flat a-Si02 substrate has not been measured di- 
rectly, but measurements of the heat transfer between 
carbon nanotubes and sapphire by Maune et al[ll| in- 
dicate that it may be of order a « 8 x 10^ W/m^K. 
This value was deduced indirectly by measuring the 
breakdown voltage of carbon nanotubes, which could be 



related to the temperature increase in the nanotubes. 
Molecular dynamics calculations [9[ for nanotubes on a- 
Si02 gives a w 3 x 10* W/m^K (here it has been as- 
sumed that the contact width between the nanotube 
and the substrate is 1/5 of the diameter of the nan- 
otube). Finally, using a so called 3 u method, Chen 
et al[l^ have measured the heat transfer coefficient a w 
2 X 10* W/m^K. 

As pointed out in Sec. 3, the model developed above 
for the heat transfer involves a single, or narrow band, of 
bending modes of the membrane or 2D-system. In order 
for this model to be valid, the coupling to the substrate 
must be so weak that the energy transfer to the substrate 
from the bending mode occur so slowly that the mode 
can be re-populated by phonon scattering processes, in 
such a way that its population is always close to what is 
expected if full thermal equilibrium would occur within 
the 2D-system. This may require high temperature in 
order for multi-phonon scattering processes to occur by 
high enough rates. We suggest that this may be one 
reason for the decrease in the heat transfer coefficient 
observed for the graphene — a-Si02 system below room 
temperature[12| . 

In Ref. [8| the temperature profile in the graphene un- 
der current was studied experimentally. The heat trans- 
fer coefficient between graphene and the a-Si02 substrate 
was determined by modeling the heat flow using the stan- 
dard heat flow equation with the heat transfer coefficient 
as the only unknown quantity. The authors found that 
using a constant (temperature independent) heat transfer 
coefficient a « 2.5 x lO"^ W/m^K the calculated temper- 
ature profiles in graphene are in good agreement with 
experiment. This a is at least 10 times smaller than ex- 
pected for perfectly fiat surfaces (see above). In an ear- 
lier paper [l3J we have studied the heat transfer between 
graphene and a-Si02. In that study we assumed that be- 
cause of surface roughness the graphene only makes par- 
tial contact with the Si02 substrate, which will reduce 
the heat transfer coefficient as compared to the perfect 
contact case. We showed how observed heat transfer can 
be explained as resulting from the spreading resistance, 
rather than the heat resistance from the area of real con- 
tact. 

It has recently been suggested jl4| that the heat trans- 
fer between graphene and a-Si02 may involve photon 
tunneling 'l5J. That is, coupling via the electromagnetic 
field between electron-hole pair excitations in graphene 
and optical phonons in a-Si02. However, our calcula- 
tions indicate that for graphene adsorbed on a-Si02 the 
field coupling gives a negligible contribution to the heat 
transfer [13|. 

To summarize, we have studied the heat transfer be- 
tween weakly coupled systems with fiat interface. We 
have presented simple analytical results which can be 



used to estimate the heat transfer coefficient. De- 
tailed results was presented for the heat transfer be- 
tween a membrane and a semi-infinite solid. For this 
case the energy transfer is dominated by energy exchange 
between a bending vibrational mode of the graphene, 
and a Rayleigh surface phonon mode of the substrate. 
This model assumes implicitly that, due to lattice non- 
linearity (and defects) , there exist phonon scattering pro- 
cesses which rapidly transfer energy to the bending mode 
involved in the heat exchange with the substrate. This 
may require high temperature in order for multi-phonon 
scattering processes to occur at high enough rate. As 
an application we have considerd the heat transfer be- 
tween graphene and amorphous Si02. The calculated 
value of the heat transfer coefficient was found to be in 
good agreement with the value deduced from the exper- 
imental data. 
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